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SUMMARY . 


In the present. thij;d .par.t' of an investigation dealing 
with the inward bulge type “buckling of mpno co que cylinders 
the theory of buckling In pure -'bending i s J?ef InecL through 
the consideration of the shear BtVaTtt^ensrgy i^ored fn the J 
sheet covering, the aesuAptiori T ox“A-. revised, deflecfed shape 
of the rings at buckling, and the use of summations in place 
of some of the integrations of the preceding reports. The 
calculation of the shear strain energy 1 b based oh ah esti- 
mated variation of the shear rigidity, of th¥ buokled panels 
Of sheet with increasing compressive stress, which makes it 
possible to dispense- with the empiricism involved in. the 
procedure of the preceding papers ■ according to which the 
value of the parameter n had to - be determined .from data _ 
obtained in general Instability., tests. The. revised theory 
is applied to all the sheet-covered specimens tested in pure 
bending In the Guggenheim Aeronautics Laboratory Of the 
California Institute of- Technology 'and at the Polytechnic 
Institute of Brooklyn, The agreement between theory and ex- 
periment is found to be good, ■ r. 

INTRODUCTION " 


“In the present part III of a series, of .investigations 
dealing with the inward bulge .type buokldjig- - of ■monq'cdq-ue 
.cylinders an attempt Is made to develop further 6he^ theory 
of .general instability in purs bending,. The. basic .Ideas of 
the theory were first published .in ■ 1938 in reference A._ The 
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monoc-oque was assumed to consist - of a circular sheet metal 
cylinder relnforoed with evenly spaced rings and stringers, 
one of the stringers being located at the bottom (and most 
highly compressed) fiber of the cylinder. Ihe buckling load 
was calculated with the aid of the Sayleigh-Ritz-Timoshenko 
method. In the axial, direction the shape of the inward, 
bulge that develops at buckling was assumed to be sinusoidal. 
In a transverse section the distortions were assumed to fol- 
low the general pattern shown in figure la. Ihe dotted line 
in figure lb presents mope accurately, and on an enlarged 
scale, the deflections of a ring according to the formula 
suggested in reference 1, This formula satisfies the re- 
quirements of inextensional deformations- and of symmetry and 
represents a distorted shape in which the deflection vanishoa 
and the tangent to the deflected shape is perpendicular to 
the radius of the circular- cylinder when <p <■- <p 0 . 

In the caloulati ons the shear strain, energy stored in 
the. .sjheet. covering was negleotod. Ihe buckling load P Qr ,. 

defined' as the: force acting- upon the most, highly compressed 
stringer .at .the time when the monoooque buckles, was obtained 
as a fun obi on ^of two parameters. Ope of these was denoted 
n and defined as n/q> 0 . It char ao.t crises ? the length of the 

pe.rimet'er involved - in' buckling as shown' in figure la. Ihe 
otharr was the number in 'of; rings ‘'involved in the inward- 
bulge. - The buckling. load was minimized by differentiating 
it with respect ;td; n And setting the 'differential coeffi- 
cient equal to zero.. This equation, together with the orig- 
inal expression .for' the buckling' load - , "permitted the drawing 
of a family of curves representing the equivalent length - 
factor X, defined by the equation 

P or - » TT a SI 8l . r /V ? L» a '- . .. (l) 

as a function of the buckling index A, defined as 

A « <r 4 /L a 3 d)<BI> 0t;p /(!Sl) r (2) 

s 

Ihe parameter of the family of curves was tho, number m of 
rings involved in tho inward bulgo;, Prom the diagram tl*e 
minimum critical load of any circular mCnocoque cylinder and 
the corresponding value of m could be easily determined by 
seleoting the maximum value of X - that is, the value on 
the envelope of the family of curves - corresponding to the 
buckling index A of the monocoque cylinder. 
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Unfortunately a comparison of the theory of reference 1 
with the results of the tes'te carried out in the Guggenheim 
Aeronautics Laboratory of the California Institute of tech- 
nology (see reference 2) showed that the .theoretical buckling 
load amounted only to 11 to 83 per'eent of the buckling load 
found in experiment. Ihe. reasons for-this discrepancy were 
discussed in fe ference 3. Simple ■ calculations showed that 
if the equations for the deflected shape of the rings were 
assumed to. be valid for the sheet covering also, the Shear 
strain energy stored in. the sheet greatly exceeded the strain 
energy stored in stringers .-and rings. On the other hand, in 
all the experiments performed hitherto the panels of sheet 
were in a buckled state in the vioinity of the most highly 
compressed stringer when general instability occurred. Ob- 
viously buckling of a curved panel of sheet reduces its re- 
sistance to shearing deformations ,- particularly if the radius 
of the cylinder is small and the stringer spacing large. 
Consequently , a calculation of the shear Btrain energy can- 
not correspond to reality if baser! upon an unchanged value 
of the shear modulus G- and a deflected shape of the sheet 
represented by the formula suggested for .the d’efleoted Bhape 
of the rings. ' ' - ‘ " - ' — 

Moreover, the shearing rigidity of the panels varies 
around the circumference, those close to the most highly 
compressed stringer offering little,’ those near the neutral' 
axis of the cylinder a great, deal of resistance to shearing 
deformations. Consequently even in cases when at buckling 
the shear strain energy stored in the. sheet is email as com- 
pared to the strain energy stored in stringers and rings, 
the sheet has a decisive effeot upon the length of the per- 
imeter involved in buckling, and thus greatly influenced the 
buckling load of general instability’. Because Of the great 
shearing rigidity of the non-buokled panels it is clear that 
the inward bulge is not likely to extend to the tension side 
of the monocoque oyllnder, but it is not eksy to oaioulate 
the exact circumferential length of the bulge, or the exact 
value of n which is tantamount to it., Site effect of th’e 
sheet upon n was not considered in the theory of reference 
1, and this omission resulted in an" overestimate* of ’the cir- 
cumferential length of the bulge and a’ 'consequent under- 
estimate of the buckling load. ’ ' v , 

In reference, 3 a .revised theory was' developed, based on 
the fundamental assumption that the 'buckling ■ load can be cal- 
■ culated accurately enough .for .practical’ purposes from the 
strain energy balance in stringers and rings without con- 
sidering the shear strain energy- bt ored in -the sheet covering. 
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provided, the value of the parameter n ie known from general 
instability tests. This theory yielded the following formula 
for the buckling load: 



The values of the symbols in this formula representing geo- 
metric and mechanical properties of the monoooque cylinder 
pertaining to each test specimen of the GALOIT test series 
of reference 2, and the critical load ■ observed in experiment 
were substituted into equation (3) and the equation was 
solved -for n. The values obtained ranged from 3.12 to 6.59, 
The variation was .attributed, to the effect of the parameters 
d/r and - . characterising the shearing rigidity of 

the panels, and the. values of n Were plotted against them 
in figufre i.0 of reference 3; The faired-in curves of this ' 
figure show an average deviation. from the . experimental points 
amounting to 4.4 percent. It was suggested that this figure 
be used in conjunction with equation (3) for the' prediction 
of the buokling_l oad of monoooque cylinders* 

In the present report the development , of .the theory is 
oar.rf.ed one step further. . The main improvement is t-he con- ' 
siderotlon of the shear ■ strain energy stored .in the sheet 
o over ins .permit ting 1 tho elimination of the empiricism In- . 
volved in the determination of the value of ' .n V from tho re- 
sults . of general instability tests, -the oal^o'ula’t.ion of the 
shear strain energy .was made possible by the' assumption of • 
an expression conneoting the effective shear modulus .9’ 
of a panel with the normal: 'stress ’prevailing in the stringors 
adjaejent to the panel. This effective shear, modulus, multi- ■ * 
plied by ono-half the .square" of the average shea.r strain in 
the panel, as calculated from the displacements of the four, 
corners of the panels was. taken to - represent, the shear strain 
energy s-tored. in' the panel. - A eScdnd improvement was [ the 
choice. of a new equation for the distorted shape of the ring^; 
( see fig, ..lb, ) It^ satisfies the- additional Requirement t hs*t . 
•the'fadius of curvature of the deflected .ring^-be equal • to 
the original radius of the cylinder when - cp ° qp 0 » 

/ •* • ‘ * " I - 

, ,, ■, * , j * 

In the . development of the, theory. e.Qme: of the- original 
integrations' were replaced by. summations. Consideration was 
given to the possibility of deviations of .the .normal stress 
distribution from the commonly assumed linear- law. The- re-* 
suit's .of tho- theoretical calculations are presented in the 



NAOA IS. No. 968 


6 


form of tables the use of which permits the evaluation, of the 
minimum critical load for m.onocoque cylinders having the 
same number of stringers as the test specimens of the CA'lCiT 
and the P.I.B. series. Tables also could be sot up e&ffily 
for cylinders having any other number of stringers. It is 
considered advisable, however, tp delay tho presentation of 
such additional tables until the time when more is known 
about the variation of the effective shear rigidity of ~tfie 
panels and the actual shape of deflections at the moment 
when buckling begins, and when more data are available to 
substantiate the conclusions of the present theory of the 
inward bulge. type general instability. 

.The present theory was applied to the pure bending test 
specimens of the GA.LCII series and to those of the P.I.B. 
tests described in reference, 4. G-ood agreement was found 

between theory and experiment. 

» 

This Investigation, conducted at the Polytechnic Insti- 
tute of Brooklyn, was sponsored by and conducted with the 
financial assistance of the National. Advisory Committ.ee for 
Aeronautics. 


SYMBOLS 


parameters characterizing effective shear rig- 
idity of a panel 

cross-sectional area of a stringer 

parameters involved in describing variation of 
direct strain ■ .... ■_ - 

kth coefficient of a Fourier series 

stringer spaoing measured along circumference 

Young’s modulus for material of a ting 

Young’s modulus for material of a stringer 

function of number of ring fields involved in 
failure 

function of parameter n arising in connec- 
tion with ring strain energy 


a o * a i 
A 

*o.*i 

c k 

d 

3 r 

■®str 
f (m+1) 

f r (n) 
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^ i sh 
fa flh(s,n) 
Fgh ( a , n )• 

f str r ^ » n 
f str t ^ a * b- 
2*8 t r ^ R » n ^ 
£y(s ,n) , 

g 

O' 

h 

i 

*r 


*8tr o 
*str r 
I atr t 


L functions of parameter a a and tt appearing 
f in sheas' strain energy term 


a 0 f i eh^ 8 * 11 ^ ** ft ifa sb^®*h) 


functions of parameters a and n appearing 
j j in stringer strain energy, term 


*str r ^str r^®» n ) + *str t ^etr t^ a » n .) 


function of parameters a and n appearing 
in work term 

ordinary shear modulus of a panel 

effective shear ' modulus, of a panel 

di stance "’of centroid of a stringer from sheet 
covering of oylinder 

index used in summations or products 

moment of inertia of ring oross soation and 

its effective width of sheot for bending in 
radial direction 

moment of inertia of ring oross section without 
its affective width' ‘of sheet for "bending In 
radial direction 

moment of inertial of stringer cross section 
without its effective Bheet for bending. in 
radial dir eet'l on "• *' ' * ' “ ■’ 

moment’ of inertia of stringer .oross section and 
its effective width for bending in radial 
direction" _ 


moment of inertia of stringer cross section and 
its effective width for bending in tangen- 
tial dlre-ctton 



NAOA TK No. 968 


r 


Istr t o moment of inertia of stringer cross section without 
its .effective width for handing in tangential 
direction 


* 


•« 


* 


i 

k, z, 

L 

(m+1 ) 
M« 


n 


P 


P 


or 


index used in summations or products 

integer constants 

total length of wave in axial direction ' 

distance between adjacent rings measured’ axi^tly 
along cylinder _.i_ 

number of ring fields inaluded in length L 

function of (m+l) associated with shear Btrain 

energy term - - 

a parameter cbaraoterizing shape of ring deflection 
which is equal to ratio Oif total circumference 
to that inv-olved in buckling, or n « Tr/cp 0 

a parameter which affectB the direct Btrain dr 
shear rigidity distributions 

load carried by most highly compressed stringer 
and its effective width of sheet at buckling 


pi 


load carried by any stringer and ltB effective 
width of sheet at buckling.. 


r 

E 


s 


■S 


*“t 

u 

U 


radius of circular cylinder ■ 
polar coordinate 

number of stringer fields Included in bulge 
total number of stringers in structure 
thlokness of sheet oovering 

integer constant . . 

strain energy 


U 


r 


strain energy stored in rings 
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tTgjj strain energy stored ir. sheet covering 

« 

^str strain energy stored in stringers 

V S = (4/3) rrr 5 , volume of a sphere of radius r 

2w 0 ordinal effective width associated with most highly 
“* compressed stringer 

2w‘ * modified effective width associated with a stringer - 

2w* average modified effective width of all Btringors 

w r radial displacement of a poi-nt on a ring or a stringer 

wt tangential displacement of a point on a ring or a 

stringer 

V virtual work done "by external load 

x coordinate measiiring distance along axis of cylinder 

t 

y * deflection of a straight beam due to heading 

16a maximum deflection of most highly compressed stringer 4 

or most highly deflected ring 

16P maximum deflootion of any ring 

y average angle of twist of a panel . 

6 shift of neutral axis measured from horizontal geo- , 

metric axis in terms of radius 

AL ohange in distanoe between x- » 0 and x = L due 

to distortions 

€ max strain in most highly compressed stringer at failure 

* *k 

A equivalent length factor k 

A buckling index . , 

p radius of curvature at any point of a ring 

cp angular measure 
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q> 0 ' angle defining end of. cirqumf efehtial ware 

dty change in direction of tangent 


DEVELOPMENT 03* THE HUT IS ED THEORY " ' " 

Revised Deflection Pattern of the Rings 

She following seven requirements for ( tlie distorted 
shape of tho rings were stated in reference 1, though in a 
slightly different form! 

' 1 , — ■ 

(1) Distortions must he inextensional. This require- 

ment can he formulated as 

w r = - dw^/dcp (4a) 

where w r ■ and w^. are the radial and the tan- 
gential displacements, respectively, of points 
lying on rings or stringers. the proof of this 
statement is now. given w^h the add' of figure 3. 
The part of_ the el ongatl on of thcTarc element 
ds = rdcp that is’ due to a radial displacement 
of .tho element is’ 

[w r + (l/2) (dw r /dcp)dqp]dcp 

The part that Is caused hy the tangential dis- 
placements Is - ‘ 

w t + (dw t /dcp)dcp - w t 

The requirement that the length of tho arc ele- 
ment remain unchanged can, therefore, he ex- 
pressed as - - - 

« . * 

w r dep + (dw t /dcp)dcp a 0 (a) 

if a term containing th_e square -of dqj is 
omitted as a .second order small quantity. Prom 
this- expression equation (4a) follows immediately 
The fact that curved "bars usually distort Inex— ' 
tonsionally . was- first noted hy Lord Rayleigh. 

(2) The maximum radial displacement should occur at the 

location of the most highly compressed stringer! 

dw r /dcp = 0 when cp = 0 (4h) 



NAOA. TN No, 968 


The t angjani i a'l displacement must vanish at the lo- 
cation of the most highly compressed stringer: 

wt = 0 when qp «■— 0 (4o) 

The radial displacement' must vanish - when qp = cp 0 t 

w r = 0 when' cp = cp 0 (4d) 

The tangential displacement must vanish when ■ 

<P - <P 0 * . 

i 

w^. » 0 when qp = cp 0 (4s) 

The tangent to the deflection- curve* must be perpen- 
dicular to the radius of the nondistorted cylin- 
der at the -end of the bulge. This requirement 
may be expressed as 

dv r /dqp *- 0 . when- qp «= qp. 0 . (4f ) 

The proof of thiB statement is now given with 
. the aid of. figure 2, The change irT the direc- 
tion of the infinitesimal arc element io meas- 
ured by dqpr - It fallows from the -figure that 

tan(d\J/) 1 = • [(dw r /djp;)d’q)]y [> +. w r + (dw t /dqp) ]dqp 

In -the limit .the tangent can' be replaced by the 
angle. Because of equation (a) the denominator 
reduo^js to r. Consequently 

d\}/ = (l/r ) (dw r /dqp) (b) 

The radius of curvature of the deflected shape at 

the end of the bulge must bo equal to the origir=*- 
nal radius of the cylinder. 

Since the expression fur tho curvature (l/p) is 
in polar coordinates 

(1/p) = [R a + C (dR/ dqp) 3 - R (d 3 R/dq> 3 ) 3/ [R a +■ CdH/d<p) 3 ]% 

and since the radial distance R of any point of 
the distorted median line of the ring fro* the 
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origin. Is the sum of the original radius r and 
the radial displacement Wr» the curvature 6 'afi ' • 
he written as 

(l/p) «= (l/r a )[r-w r - (d a w r /dcp a )3 (c) 

after substitutions are made and powers and prod- 
ucts of v r and its derivatives are omitted as 
second order small quantities. The requirement 
for an unchanged curvature is "■ 

r - w r - (d w r /dep ) = 0 

This requirement Is equivalent to the following 
condition, if requir einent (4d) is fulfilled: 

d w r /dqp a = 0 when cp *= cp 0 (4g) 

The expression assumed for w r . in the preceding reports 
satisfied all hut the last fcequir emefctV- Nonfulfillment "dY ' 
equation (4g) is believed to-haye had little eTfect upon the 
buckling load. In the present Calculations a" hew expression 
for the deflected shape is used which satisfies all the sevSn 
condi-tions. It is found- that the shear strain energy corre- 
sponding to the new expression -is materially Smaller than 
that corresponding to the original expression. The new ex- 
pression is as follows: > - 

w r = - j3 [2 cos(ncp/2) + 9 ooe(3ncp/2) + 5 cos (Bncp/s) 3 (d) 


when 

and 

when 

In equation (5) 


’“(tt/b) < cp <. (Tr/n) 
w r = 0 

. jqp,j > (tt/u) .... 

3 *= (a/ 2 )[l - cos (2tt3:/Ii) 3 



(db) 


where a is an indeterminate coefficient proportional to 
the maximum deflection of the cylinder. It may be seen that 
g represents the same function, it represented in reference 
1. Consequently the shape of the inward bulge is unchanged 
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in axial sections of the monocoque cylinder. Because of the 
requir ament given in: equation (4a) the tangential displace- 
jaent w^. can he written as 

w^. = (0/n)[4 sin(ncp/3) * 6 sin(3nqp/3) + 3 sin(6nqp/2)3 (6) 

-(-rr/n) <_ <p < (u/n) 

= 0 (6a) 

Jcpj > (TT/n) 

It is easy to prove that equations (6) and (6) satisfy 
all the conditions presented in equations (4), The old and 
the new deflected shapes of the rings are compared in figure 
lh for the case when n = 3. 

Obviously the choice of equation (6) is arbitrary. 

Other different equations can be found which satisfy~'the 
seven requirements. It is shown in appendix I how expres- 
sions can be determined which, while fulfilling the condi- 
tions set forth in equations (4), satisfy additional require- 
ments and correspond to -further reduced values of the shear 
strain energy. These - expressions represent deflection Curves 
that hug the original circular shape more closely in the 
neighborhood of cp 'tt <p Q -, as - may be seen fyom the curves 
given in figure 3, 


when 

and 

when 



Strain Energy Stored in the Bings 

The strain energy stored in a ring oan be calculated 
from the equation 

(l/2) (EI) r [(l/p) - (l/r)] 8 r d<p (7") ■ 

where * 

(El) r bending rigidity of ring for bending in the plane of 
ring 

r original radius of cylinder 

p radius of curvature at any point after distortions 




i 


i 
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The moment- of Inertia of the ring is calculated for the ring 
section proper augmented with its effective she'et. The width 
of the effective sheet is. taken equal to the width of the 
ring in agreemeh't with the procedure advocated in references 
1, 3, and 5. With the aid of equation. (c) the strain energy 
stored in all the rings can he written in the .form! 


+Tr/n 


U r = V (1/3) (ll) r / (l/r 4 )[w r + (d a w r /&cp a ) ] a r dqj (8) 

U 'W( a 


where the summation must include all the rings involved in 1 
wave length,- Upon substitution, of the value of w r as given 
in equation (5), and with the aid of the transformation 

x = JL/(m+l) 

valid at the location of any ring, integration gives _ 

m ' " . ..... 

U r a (aa/4)TTf r Cn) [(ai) r /2r 3 3 ^ ^1 - cos [3ttJ/ (m+1 )]J* (9) 

3=i ’ ’ ' ‘ 

where ' ' r ' ' ~ : 


TT 


/n 


'* fr(w) = Cl/W) J* [W r + (* 3 w r /bep 3 )] 3 dcp 

-tt/u 


(9 a) 


In general, if w r is given hy the expression 


w r ~ ” ^ v c i oos kj_ncp (2k^ a positive integer) (a) 

i-i ,a ,s, . . . 

then . equatl on (9a) • beoomeB 

fy(n) * (l/n) 

*’ i 


. , (B| 
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With the values of the coefficients c given in equation 
(6) equation: (B) eimplifies to 

f r (n) =[1386.875 ni - 679 n + (llO/n)] ClO) 

The summation indicated in equation (9) was evaluated 
in the appendix of reference 3. It wbb found to he equal to 

(3/2) (m+l ) when (mil) > 2 (lla) 

1 when (m+l) = 2 (llh) 

Thus the total strain energy stored in all the rings within, 
one inward bulge can he written in the following form: 

U r =a a ( 3 /I 6 ) [ (m+l )/r a ] (-rr/n) [1386.875 n 4 -679 n®+110 ]'(SI ) r (12) 

when (m+l) > 2 

By replacing (3^2) (m+l) hy unity in the right-hand member of 
the equation, the. expression for the strain energy valid 
when (ra+l) => 2 is obtained. 

It is seen that the total strain energy stored in the 
rings is proportional to the. number of r-lng fields (m+l) in- 
volved in one bulgo. Moroovor, tho strain energy Increases 
very rapidly with increasing n^ approximately as the cubo 
of n - that is, it takes very much strain energy to distort 
the rings into bulges of small circumferential length. Tor 
numerical caloulati-ons it is convenient to rewrite equation 
(12) in the form 


u r = a a (3/l6)'[Cm+l)/r 3 ]TT(El) r f r (n) 

when (m+l ) 

.>' 2 

(13) 


U r = a 8 [l/ (8r 3 ) Jtt( 2I ) r f r (n ) 

when (m+l) 

= 2 

(13a) 


where f r (n) is the function defined i 

u equation 

( 10 ). 

7al- 

1 

ues of this function are tabulated in 

table i. . 

— 




Strain Bnergy Stored in the Stringers 

In accordance with the assumptions made regarding deflec- 
tions, all the stringers band both radially and tangentially 
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(except the most highly compressed- stringer which. bends,. only 
radially), The strain energy oah be calculated. from equation 
(7). It ie shown in . Strength . of . Materials that for an origr 
inally straight 'beam ‘ -- -- . 

(l/p) - (l/r ) « d a y/dx a . , 7 ( 4 ) 


if .y designates the displacement perpendioular to the 
x-direction (x heing measured parallel to the axis of th6_^ 
"beam). As w r and w^' denote two perpendioular components 

of the displacements of' the stringer, they have to he UB'ed 
in the place of y. . Thus the totals train energy stored in 
the stringers can.be written in. the form t~'~\ 


u atr = ft/S^str J*. V^^etr, t + (a 3 w r /3X P ) 8 I atr r] ( 


dx (lb) 


. 0 . . . . . .. . .. -V. T 

where the summati on :must include - all the., str ingers • It is, - 
assumed that. I Btr ., t . and I gtr r ,‘ - are. .principal moments of 

Inertia. • After substitution of the expressions; given in ~ 
equations (5)'and (6) into equation ('14) integration glvoa 


u str = (a-V/L 3 ):ili str 'Y '{k C08 ^ q ¥ , / 2 ^ 4 '9 oosC3np/ 2 ) +' 5 cos(5rqV 2 >3 a I 8 tr r 

» ““ ' \ T # - * r - * *-*- is #1 f* t¥tt^ 

+ (l/n 3 ) {h Bln(ncp/a) + 6 sin(3ncp/ 2 ) + 2 sln(5rcp/ 2 )3 I B tr tj ^5) 


The summation has to he extended around., the part of the per- 
imeter Involved in buckling. In this region l a tr r ail "^ 

^str t vary since the effective width of. sheet added to the 
cross section of the stringer proper varies because of the 
changes in the value of the normal stress. 'The analytical 
expressions for the moments of inertia are’very cumbersome 
particularly since the effect of the curvature of the sheet 
covering must also be- taken into ao.count. Moreover, the so- 
called reduced .effective width ratlidr than ’the ordinary ef- 
feotive width muat .be used f or • tho .■calculation 0 ? ^6he moments 
of inertia, as was • shown in referonoe 1. ' Assuming Karguerre’s 
cube-root formula I^Q. be .valid for the effective width 2w of 
the equivalent fl.’.t sheet, the reduoed effective width , 3w* 

— (s/3)2w in the region where the defleotions are in the In- 
ward. direction. In the region- whore the cylinder bulges 



16 


SAGA SIT So. 968 


"out” (see fig. la), the "reduced^' • effective width is 4/3 
of the ordinary effective width and should "be denoted prop** 
erly as " "modified" effective width. -This term obviously may 
refer to both the deoreased and increased values. 

The complexity of the problem, can. be reduced materially 
through the assumption of average moments of inertia ^etr r 

and ^atr t corresponding to the average value, of the modi- 
fied effective width. The rule according to which this av- 
erage must be choson cannot be determined, howevor, without 
carrying ottt the; summations of equation (15) with all mathe- 
matical rigor,' On the other hand, the critioal strain in in— 
.ward bulge typo buckling is not too sensitive to changes in 
the rule usod for doteraining the average value of the modi- 
fied offoctivo width, as may bo soen from tho numerical work 
of tho fallowing chapter, and ovon les.s so to changes in tho 
choice of the value of the effective width of the most highly 
compressed stringer, as was shown in reference 6, Moreover, 
it ehduld be remembered that the conception, of the offeotive 
width itself is approximate and in particular the applica- 
tion of the conception to curved panels is laoking in ana- 
lytical foundation. Consequently great rigor is not called 
for in the calculation of the effective width, 

Tor these reasons the theory of this report has been 
developed on. the assumption of a oonstant average value for 
tho moment of inertia. In all tho numerical work to follow 
tho average momonts of inertia wore caloulatod in correspond- 
ence with the- ostimate that the value of the -averago modifiod 
offoctivo width lies botwoon the values of tho ordinary and 
tho roduood offootivo widths at tho location of tho most 
highly compr-ossod stringor, Tno principles govorning tho 
choioe of the value wi-thin these limits are discussed, in con- 
nection with details of the numerical work. 

The assumption of oonstant moments of Inertia permits a 
vigorous calculation of the summation indicated in equation 
(15), The work involved oan be slightly simplified through 
the introduction of the symbols S for the total number of 
stringers (or stringer fields) in tho monocoque cylinder^ and 
s for the total number of stringer fields affected by-' the 
distortions. With this notation the angular spacing of the 
evenly spaced stringers is Sh/S radians, Moreover, from 
the definition of n 


s ** S/n 


(16) 



TJACA TN Ho. 968 


17 


With the new notation equation (15) can f® written in 
the form . * • ' - 


U g ^r (L 3 /a a Tr 4 )(l/® s tr) = S’str* 8 ’ 11 ) = X str r f str • 

• •. Xs/aJ-x 

- + Wt^tr t (s * n) ' 256 Z str r + 2 


[2 cos (rri / a ) 


1 = 1 


Tl7) 


2 

+ 9 cos(3iTi/s) + 5 ' cos (5Tfi/ S ) ] ^gtr r 

( S/a )— i _ .... _ V - 

• +■ (-2/n S ) ^ [4 ain(TTi/S) + 6 iln(Snl/S) 

i = i. 

•+ 2 sin(5TTi/S)J Ig-fcr t • . 

The term 256 I str r represents the strain energy stored in 

'the mopt highly compr.es sed s tr inge^ . * s tUingWB. 

resent the strain energy stored in all the other 
These summations are carried out in appendix II. They give 

V 8 tr (s.n) ='.55 8 I 8tr ' r > 38(s/n a ')l tftV . t when a >4 (18a) 

r str (s,n) = 400 I Btr r + (64/n 2 )l si?r t 
T’ st . r (s ,n) = 256 Igtr r' 


when s=4 Tl80l 

when s = 2 (18c) 


The 'functions £ str r (s,n) and f str t^ s » n ) are Tab- 
ulated in table. 1. With the aid of the function ^tr^> n 
defined in equations (17) and ( 18 > the^.tr*ln ejeW^tored. 
in all the stringers can be given in the concise 

(19) 


U str '.= (a, 3 Tr 4 /L 3 )B str 5 , str (s ,nT ■. 
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It is interesting to note that calculation of the strain 
energy stored in the stringers through integration, based on 
the assumption of continuously distributed stringers, gives 
the same expression as that obtained by the preceding summa- 
tion, provided the number s • of stringer fields involved in 
buckling is eq.ual to or greater than six. It also may be 
seen that as the - ring strain energy was found to be directly 
proportional to the number of ring fields, so is the 
stringer strain energy directly proportional to the number 
of stringer fields involved in distortions (except when the 
number of stringer fields involved is less than six). 


Shear Strain Energy Stored, in the .Sheet Covering 


It was mentioned in the Introduction that the main dif- 
ference between the approach of the preceding reports and 
tha_t of the present paper is the consideration' in the present 
paper of the shear -strain .energy . stored in the sheet cover- 
ing. The calculation of this shear strain energy was made 
possible through the use. of an average valud Of the shear 
strain and of an average value of the effective shear modu- 
lus for oach panel. Figure 4 shows a panel before and after 
■distortions. With the notation of the figure * fe.be .average 
angle of shear of the panel may be expressed as 


. Y 


* ■<1/3L 1 )| (w t V i<k .,.-+ (w. t ) (cpi+ ( 2Tr / a; ). Xl ) 
:: ~-J^ w t\p^,x fe + fat ^.(cpj+CaTt/s •) ,x 3 ■);. f 


(30) 


The . t otal .. shear strain energy stored in all the panels is 


U 


sh 




(1/3) v Or* Lj td 


( 21 ) 


.where . .the summation -has: tp be . carried out, over all the panels 
involved, in the bulge. In this .equation stands for the 

effective shear modulus of the panel. When the panel is in 
the non-buckled state, ' 0} ' is; -e.q.u© - !.;. &o ? : • If, hovovop, 

the coaprosaive’ stress in the panel. 1 -exceeds.; the critical 
stress, it is reasonable to assume that the panel will offer 
less resistance to. shearing di.s.t orti ons .. The decrease in 
the shoar rigidity is likely to dopond both on tho ratio of 
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the average compressive stress -in the panel to- the -critical 
stress of the panel, and on the ratio of the stringer spac- 
ing to the radius of curvature Of the cylinder, A theoret- 
ical evaluation of this dependence appears ( to "be very diffi- 
cult, An experimental est Ahlishment • of the connections 
should he easier, hut no publication dealing with the sub- 
ject is known to the authors. Because of this lack of in- 
formation it was necessary to assume arbitrarily an expres- 
sion for the effective shear modulus. It was decided to use 
the equation 


. G 1 


. G 


a 0 - . aj, (cos cp) 



( 22 ) 


If In this expression a 0 , ax, and p‘""ar6 equal to unity, 

the value of the effective shear modulus .Is . zero at the most 
highly compressed stringer,- and ■■ G at the horizontal diam- 
eter of the cylinder. If a Q is equal to, and aj slightly 


smaller than unity, G* does not become zero even near the 
most highly compressed stringer, -In. both cases the varia- 
tion of ~G ,- ' is linear' with vertical "distance 'from the most 
highly compressed stringer, ‘This linear law can be changed 
to a parabolic one by assuming p to bo greater than unity, 
Figure 5 shows theso possibilities. 

It is reasonable to assume that ; other 'Chingpi^elhg ~ 
equal, the value of G* will be small in panels upon which 
large compressive forces are acting, and large in panels 
upon which small compressive forcee are acting, Thus the 
linear law for G* ^ust discussed may' well represent the 
case of a cylinder in Which the normal stre"ss distribution 
is linear. The assumption of a parabolic normal strain dis- 
tribution shown by the dotted line in figure 6 appears to bo 
in, reasonably good agreement with the experimental curve 
shown by the full line in figure 6, if the region alone is Q 
considered in which the bulge, develops at buckling (g> 0 < 90 ) . 

The curve represents the normal .strain" distribution ix^ speci- 
men 159 of the GALCIT series and is taken" from reference 6. 
Consequently -aquation (22) with p = 2 might well represent 
the variation of G l around the circumference of a moftdeoqus 
.cylinder in which the distribution of the__n-ormal- stress, is 
not linear, and the, peutrai axis is materially’ shifted from 
the center of, the cylinder. 


In any case, it is hoped that ' equation" (22) c.anbe well 
enough adapted to experimental curves which might be avail- 
able in the future* : Calculating from equation (22) the - 
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panel “‘•ear nodulu. for the. middle of a. , 

. •' * °f 0 -- ..[to, ^.li±2lllp\ (23a) 

: ° f a ^ a - ti °*e (6), (20 T, and (22a) into" - 

equation. (31) there results after s one . s imp! if icat i on 

m . (#'/s )—l J 

T I'"' .I 

: '.j=o .. i=o . . r J 

J * 3 [sin, S2i + "tin Milii ] V 


sh = 


cc 3 0t d 
32L- 


. + 6 


si tf&Li. + gin g-Tt Ci+1 ) 

•'■ • ts / •' ■: ' s 


{ a 0- -®1 f 0 03 ] v p, 1 r coa - + CQ 

1 " ■ -L • : ... .. .S J JL m+1- ■ CQa b+T J 


8 ' J J 

(23) 


SslT 


The summation over the ring fields i 
^,1* . With the notation . , 


s carried out in appen- 


j = o 


grr(.j'-fl) 


• . m+1 


COB 

J ' . ' ja+1' 


the result is 


f (m+1) 
m+1, 


(m+1) 


(i. 


= 8 


2tt \ 

003 m+T/ • when -Ja+1. >. 2 
• ■ when?- m+1' =' a 


(34a) 

(34h) 


!ithout°iimit f appr0 ^ cixe8 , . ^ *hen ' n+1 ■- increases 

without limit. It is graphed i n , figure; ?, - <• -V r . ■ - 

flrst ^-p :aPP ^ i,X 11 the product, of th ; e. expreesion in the 

otll ° f equatl0n X23;) .b> ‘ ,a 0 .was summed up 

over all the -stringer fields . I), sing the .notation ^i x th („ ta 
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’.for this sum. Write 


ff'4* 




a 0 f a. s } 1 (s ,n)=(gh 0 /n )s 56+4 [4 cos, ( tt/s )+9 cog (3tt/ s ) + c.os (-Stt/ s ) ] 


■ 1; 


_t - 


provided • s > 4 


~ C25) 


It was not found practicable to sum up t.h-e product of the 
same expression "by the second term'in the second brace. It 
was easier to calculate it numerically, for. the different •• 
numbers of stringers contained in jfche mono.Q, 04 ue. djfTTn'der and 
for the different numbers of stringers Involved in buokling. 

-The numerical values of these functions . a-ye-, tabulated in- 
table 1^ Using the notation’ -a x f g . ( s-* n) f.or the second 

function, write ' ' - 4V ri *■ 


U 


sh ” 


= a 


G-td 

321, 


f (m+l) 
m+1 . 


t 


a Q ’f i sh(s,n) - a^ashCs, n) 


(26) 


liquation (26) can be written In a more concise form by using 
the function I’ g j l (s,n) dofinod as ^ • 


P sh( s ’* n .) * * 0 f i ghCs/nJ- - a x f 3 g h Cs,n) ~(26'a) 


Work Done by the External Forces 

The loading of the monocoque cylinder consists of a 
pure bending moment which Is resisted by the stringer plus 
effective widthof sheet' Qomb inat i one . Consequently the 
work Qf the external forces is the sum of the products of 
the normal force acting upon each stringer plus effective 
width combination by the change in the length of the stringor 
corresponding to the assumed defloctod shape at 'buckling. 

In the proceding reports on the inward bulge type of general 
instability the stringers were assumed to be continuously 
distributed around the perimeter, and the work was calcu- 
lated by integration. The accuracy of the theory is' now in- 
creased by treating the stringers individually, which proce- 
dure, of course, necessitates replacing the integration by 
summation. "• r 

The normal force acting upon a stringer pluf, 

f eqtiye sheet 
equation 


combination is assumed to b$^ girenr by the 


P* K * b-i 


a» 

,b 


v 


(37) 


• +• 
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where -^cr the normal force acting upon the. moat highly 

compressed stringer plug effective sheet combination when 
buckling occurs. The form Of equation (27) is similar- to 
that of equation (22). If the normal force is distributed 
linearly and the symbol S is used to denote the ratio of 
tho shift of the neutral axis to the radius. 

^o = . ®/.£ 1 + 5 ) and. bj = 1/ (l + s) t provided p «« 1 . 

A nonlinear .normal force distribution can bo reprosontod by 
tho sa/mo formula if p > 1, Pigure 6 sorvos to- illustrate 
some of these ! statements. The shortening, .due to bending of 
tho distance between the end points corresponding jbo x"* 0 
and x = L of an originally straight bar is given by.. , 


AL 


(dy/ dx j 3 dx 



The effent of displacements in the tangential and in the ra- 
dial direction is additive* ' Con-s aquently the work done by 
the external forces can be written in the form 


v- y p»al 

l 

, /n n 

= (p 


r l/rf ^ J j^bo+bi (cos cp) ^ 

■ o 


(dw r /5x) + (d 


v t /bx). 3 'jij 


(28 


where the summation must include all the stringers involved 
in buckling. By using equations (5) and (6) . the following 
equation is obtained after integration ■’ * 


V = 


a a Tt 3 P 


4L ^ / 256 (b 0 +b a ) 
-A \ 


+ 2 J{[' ,o+ 1 ,i ( 00 * 2 i i ) p J [( 3 


COS— + 9 cos^i + 5 coa^Si.'N 
3 v J 4* / 


+ (l/n s )(4 sin^i- + 6 sin^SA- + 2 sin^^ 


]}) 


(29) 
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The summation of the terms containing b 0 ■ leads to ex- 
pressions similar to those given in equations ( 1 8 ) . . In the 
case of the terms containing bj. it is more convenient to 
carry out the summation numerically for given values of a- 
and S. With the notation 


* 

(S/2 )-l 


f w (s,n)=256+3 


2 cos^- + 9 cos + 5 oos^ 1 

6- S S 


I. {(••■¥)■[(' 

1=1 

a 

(l/n ) ^4 sin^p + 6 s i « - + ' 2 ~ s.i ^ 


¥) 


(2.9 a.) 

the work done hy the external forces may he given in the form 


W = 55 + (28/n 3 )]sh 0 +' b^W^ 8 ,n) (30) 

when, s > 4 


When’ s is equal to or smaller than- 4’, equation (30 ) must 
he modified. The modified forms can he written down without 
difficulty with the aid of the. expressions given £n equations 

(18). ' •"-= 7" 


A further simplification of the work expression is pos- 
sible through a consideration of the effect' upon the buckling 
load of different normal stress .distributions^ It was shown 
in reference 4 that deviations from the linear law caused 
but small changes in the work done by the external forces. 
Moreover, it follows indirectly from .the considerations of 
reference 5 that- the effect upon the .external work of a 
shift of the neutral axis may be similarly neglected. It is 
reasonable, therefore, to simplify the calculations of the 
external work through the assumptions; of linearity and zero’ 
shift - that is, to assume that ■ 


6 = 0 . b 0 = O’ “ bj = 1. ' p = 1 ( 31 ) 

On the other hand, any similar simplifying assumptions in. 
the expression for the reduced shear modulus (equation (22)) 
lead to considerable changes in the shear strain energy terms. 
Although it is admitted that there is some. ’Connection between 
equations (22) and (27), it is not considered ’ inconsistent to 
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keep the . general law for the reduced, modulus while simplify- 
ing the normal stress expression through the ' as sumpti ons con- 
tained in equations (31 ) "because the exabt' nature of the var- 
iation of G * . with the normal stress is not known at pros- 
ent. 

If, however, oquations (31) hold true, oquation (30) may 
"be written in. tho simplifioi form .... 

W »- (a 3 TT 3 /4L)P ar [f w (a,'n)3 <30a) 

This equation is valid for any positive even value of s; s 
cannot assume att odd Value because of.;the symmetry of the 
"buckled shape Vith respect to the "bottom stringer. Numeri- 
cal values of the function fy(s,n) are listed in tahlo 1. 

' v . . 

r .Determination of the Buckling Load 

The "buckling criterion can "be written in the form 

tf'tot - W = 0 (32) ‘ 

By equating the; ,.expr e.s si on of. the external work given in .. 
equation (30a) -to the sum of the right-hand, members- of equar 

tlons (13 ) , (1:9), and (26), the following .condition -for P cr * 

is obtainodJ 

*( a 2 rr z J 4L ) P^ f ^ ( s ; , n ) " i -_‘ " ' " ' . .' 

» a a (3rr/l6-5(m+l) C(El) r /r 3 ]f r (n) + (a s n 4 /L 3 )P 8 tr ( s , n )£ 

’ * ’ . T i " 

+ U a /3£) (Gtd/i. 1 ) <^[f (m+l)]/(m+i)^. P s h(B,h) • / '(33) 

. , . ! ’ * 

This equation. can-be simplified slightly through ' the. assump- 
tion that all the elements of .the. raonocoque are made of the 
same material. Dividing through by 3i(A + 2%t ) , whore 2w^ ■ * 

is- the total (not the reduced or modified) offectivo width 
of the most highly compressed stringer, it is possible to 
solYe for. the: strain----; ? cr ; in. the. most highly compressed * 

stringer at .buckling,. In the divisor t^ie effective, width : 
and not ■ tho reduced- .or modifi ed effectiYO- width haj? t 0 
used in' accordance with tho statements, mado in roferonce 1 
in connect ion with- this problom, -. ■; 


i 
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And 
€ cr 


With. Vs = (4/3)fTr 3 the volume of a sphere of radius 

d = 2TTr/S the perimetric distahhe between 1 stringers * 

2ji = l/-(m+l )■ . the ring spacing, the critical strain 
becomes af t er ’ s 6m : e' manipulat ion / v\ .. •• \ . 


■max “ ^^I>i7Tg)f T (n) (in+l) a i r . + C 2 tt 7 S r (s,n)/ (m+l) 

+’ 0..0.307 (tr7.S)f (m+l)F sh (e ,n) l/. [A+2w*t)f w (s,n)3 (34 ) 

« * * A ‘ . * _ % . k ^ r , 

i - 

According to_ equation (34) the critical strain € cr " is 

a ‘function of. . the .two .parameters' . (m+l) and ■These param- 

eters must assume integral values '' thab 'malee : 'the critical .■ 
strain a minimum. The critical s t‘r ain/cah^ho Al'nimized with- 
out difficulty with respect to (m+i) - if"' th6 &3 sumption . is 
made that (m+l) may assume any (integral or non-integral.) 
value. The critical strain then "becomes a: continuous . func- 
tion of- the parameter (m+l) so that the minimum condition 
can "be obtained' through’ different iationiy - Hence - 

ae max/ a ( m+1 > =' 0 2(L 1 /V S ) (m+l-)f T .(.n)l r . . . ./ ~~ 




-2.(2tt/, 1 1 ) 5’ gt? ,_(s f nl / Timely 
+ 2 ( 0. 03 07 ) (-hr/ S )<* (.m+l ) - Cl - -cos (2-rt/m+l ) 3 


’ - it ' sin (2-rr/m+l )j ( s r n ) '/VW+>1 ^ " 

* * . ■' ' J J" _ _ i 

If the symbol. M* . r is introduced by the definition 


M* =' (m+l ) 


(m+l ) 


(i - 




IT sin 


2rr 

m+l 


] 


■V 7 ” ‘ 

.(357 


(36) 


and the’ value of • the function is Calculated for di.fTer-~ 

ent values of (jn+1 ) , ‘ the . functf on will' be. found almost con- 
stant-. This may.be. seen from table 2. .'.‘.Since M’/Cm+l) 3 - 
appears in e<quat i.on (-35.) im a’ .term whicb i‘s usually small , 
the accuracy: ,cf. ;t.he r ...calcula ; tiohs : i s' hot'' impaired by ropTacIng 
M 1 by the ffl.ona'itanti.v-a'lue 60^’ - Thts ^sUbst''i''tutio_tt reduces 
oqUati on . ( 35) ..to the- simple* fourth, degree, ■erq.ua'f i on 
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(m+l) 4 =- [(2TT/L;! ) 2 I , str is ,n) ’ 

-■ .1.8(tr/S)y sh (s 1 n)]/.{;(L 1 /V s )fr(n)l r ] (3?) 

Equation (37) in conjunction with equation (34) would 

■ represent thq Complet e ' soluti on of the . problem of the inward 
bulge type buckling of a monocoque cylinder if the correot 
value of, the parameter s were known. This value, however, 
can be dot ermihod ' rol at ivoly easily with the aid of the 
trial-and-error method suggested In the following section. 

If greater accuracy is des ir ed , 1 equat i on (37) oan be 
used as a first approximation for, .(m+l). Then the correct 
value of M 1 can be taken from table. 3 and equati on (37 ) 

■ written in the modified form 

(m+l) =* [(2TT/L!) T gtr (a,n) 

- (0,0307).ii» (tr/S)P ah (s,n)3/C(L 1 /V s ) > f r (n)I r ] '(37a) 
This was done in the numerical calculations . 

Procedure for Calculating, the Critical Strain 

As stated previously, the monocoque cylinder discussed 
in this- report is assumed to have a stringer placed at the 
bottom (cp = 0), The value of s representing the stringer 
fields involved in buckling must be, therefore, an even in- 
teger. On the other hand, deformations at buokling are not 
likely to extend to the tension side of the monocoque because 
of thp great shearing rigidity of the non-buckled panels. 

The investigations reported in reference 3 even shoved that 
<p 0 * Tr/n usually varied only between 35° and 70°. Conse- 
quently the number ' qf the values s may assume is re- 
stricted sufficiently to pormit a trial-and-error procedure 
based on simplo guessing if. the number S_ of etringors con- 
tained in the monocoque is not too large. This was the case 
with all the Q-AIiCIT and P.I.B. specimens. One, two, or 
three likely values were assumed for s , the ' corresponding 
values of (m+l) were determined from equation (37), and 
the critical strain was calculated for each of the assump- 
tions by substituting the values of the proper functions of 
(m+l) int o . equat i on (34). The 1 " smallest of the critical 
strains so obtained was considered the actual critical strain 
of the monocoque cylinder. 
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When the monocoque cylinder contains"a greater Sumter 
of str ingers there . may be t_o.o many possible \values for - s 
to permit close guessing. In such cases recourse' may be had 
to the s emi ompir i cal data presented in figure. 10 of reference 
3 and in figure 18 of reference 5. ■ Thd predicted value of 
n can be taken from the. figure if the geometric and mechan- 
ical properties of ^he structure, and thb 'maximum strain . 
e cr are known. The maximum strain in turn- may be calculated 

approximately from the equation ' • • • • - 7 



'1 s'/T" ^ v I str o?r'o 

A n J L a rs , . . V 



where the subscript o signifies that the 'moment's of inertia 
of stringer and ring are calculated without consideration of 
the effective width of sheet. With an approximate value four 
n thus calculated, the minimum' critical strain" 'can' be de- 
termined by "comparing values obtained from equations (34) ^ 

and* (37h) upon substitution of integral values of s close 
to the valxie s = S/n. . v '■ . .nr - 

To reduce to a minimum the work involved in the applica- 
tion of the procedure, tables 1 and. -3 have been prepared-, 
which enable the stress analyst to obtain, directly the nu- 
merical values of most of the functions needed for a solu--~ 
tion of the problem. . . • : . " v - 


-COMPARISON WITH jfiXPlKlKBBl . * 

t • • ■ • 

Remarks . 


The present theory was. applied to ai.3-.the- sheet-covered 
specimens tested in pure bonding at GALC-I-T., specimens 26 
through 68, the data for which are given, in reference 2, and 
the two types of. specimens tested at P.I.B. , the data for 
which can be found _in reference 4, The computations are 
presented in table. 3. . The f ollowing- remarks may" be Of help . 
in understanding details of the numerical work involved in- 
the application of the theory. 

The effective width of sheet to be added to the' ring 
section was assumed to be equal- to the width of the ring 
section, as already mentioned. The effective width of sheet 
to be added to the stringer section was calculated according 
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to the procedure presented in reference 5, The average re- 
duced (or modified) effective width was assumed to be equal 
to the reduced effective width at the most highly comprossod 
stringer when tha number of stringers involved in the bulge 
was small. '!This' ' assumption certainly underestimates the 
average bending- rigidity of the stringers, 'and consequently 
leads to a conservative value of the buckling stress* Hever 
theless, the deviation from reality should not be too great 
since the strain energy of bending is small in the stringers 
far from the most highly stressed stringer, The influence 
of the stringers upon which less than the maximum oompres- 
sive. force is acting is more noticeable when there are many 
stringers involved in the bulge. In such cases the average 
reduced effective width was assumed greater, as a maximum 
equal to the total (non-reduced ) effective width of the most 
highly compressod -stringer. 

In the expression for the reduced shear modulus G- 1 
the constant a 0 in equation (3-2) was assumed to be equal 
to unity. Similarly was taken as unity whenever S 

was equal to or smaller than 24. 3Por S = 40 the value of 
a a was assumed to be between 0,80 and 1.00, The difference 
in this choice is based on the difference in the movability 
of the buckled panels. Obviously much less movability re- 
sults from the' buckling of a narrow panol in a cylinder hav- 
ing 40 stringers than from the buckling of a wide panel in 
a cylinder having only 10 stringers, provided the radius of 
the cylinder is. the same. The ratio of the maximum stress 
to the buckling stress of the ourvod panol also has an of- 
feot upon the shearing rigidity and was approximately taken 
iniro account In the choice of a* in the case of the cyl- 
inders having 40 stringers. 

As discussed earlier, the assumption of p « 2 in 
equation (22) may represent a curvilinear shear rigidity 
di stributi on, ' pog sibly corresponding to a nonlinear normal 
stress distribution in the bent monocoque cylinder, a no- 
ticeable shift in the location of the neutral axis, or a 
combination of the two. Since in the case of test specimen# 
having but 10 stringers the neutral axis must be. consider- 
ably shifted under the bending moment that causes general 
instability, p was assumed to 'be 2 for the 10-stringer 
raonocoques. In the case of tost eylindors having 20 
stringers the criti-Oad 'strain in genoral instability was 
usually calculated on the basis of both assumptions p * 1 
and p » 2, 

The calculation of the external work term was based on 
the assumptions represented by equations (31 )• 
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Application , of the pre-sent-- theoxy to_- 'specimen 28 indi- 
cated. that the. smallest ’Buckling lo.ad would -b-e obtained when 
only one ring was' Involved “in failure. .Hence, in this .,oase 
the 'strain energy stored in the ;ring .(.col. 25.) had to "be 
cal cul'at ed .according to equation (lib) instead of equation 
(11a). : "' T " =£ 


Illustrative Example 


For a better understanding., of the numerical work pre- 
sented in table. 3 details of the calculations are explained 
using specimen. 31 as an examples R'eferende'-. 2 gives S = 20, 
Lj. = 4 inches, r,= 15 . 92 ■ inches , ‘ ■ and t =.0,01 inch. The 

ring. and stringer sections are those designated as anoT 

respectively. The first three' columns ’in -taVle ' 3 con-'; 
tain data of the specimen. Column 4 give's the value of 

- 2.467 as -taken from table. 1,1. ■ k In the fifth col- 
umn the assumed, valu'e' of s must be fi-llod An.. • As a ‘first 
trial s is assumed; to- be -4, which is--just; about -the small- 
est possible choice.' ;This gives cp 0 ••= (4/20)tt = -36°, and, ■ 
as was stated earlier, from the investigati o’hs Tof reference 
3 It follows that cp 0 ' should have a value between 35° and 
7 0° ' • * '• b -t ; • . 

The (ordinary) effective width i's apprdxlmarfcely "calcu- 
lated to be 2 inches. In estimating 2w < , the average mod- 
ified ef f ectJ^v^LVidth, the choice' lies ■ between 3 wq <* 3 
inches and (^w o 0-, * which is approximately 1.5 incffesT ’'’So- r ' 

cause of the prevalence in. the strain-energy-in-stringer 
terms of the most highly compressed stringer against the 
only one additi onal-- s tringer involved in buckling on each ■ 
side of the plane of symmetry., " 2w T = 1;7 inches appears a 
reasonable suggestion,, and hence is adopted. 

* r ’ ' 

Column 9 contains . VI . ' „ which is taken from figure 
• s t r r 

7 of reference 5,,' or, else may be .calculated from equation 
(34) of reference 5; : 


T — T 

str r ~ x str r o 


1 


h +• (t/2)- [(2w) /( 24r) 


£ 


( l/ 2wt ) + (1/A) . 
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vh'efe h -i's the distance of the centroid of the stringer 
from the sheet covering of the -cylinder. Table 1 gives 
f - r (a,n) = 400 when S = 20 and s = 4, The product of 

this, value .."by ■. T' s ^ r is filled in, in column ’10. . In o.olunm 
11 I str t is listed after calculation from equati on , (3 6 ) 
of reference 5: 

I s tr t = Igtr t‘ o + (1/12) (2w) S t 

where: ■ I gl ^ r ^ 0 is' the moment of inertia of the stringer . 

Without the effective width -of sheet for-bending in the tan- 
gential ■ direction. Multiplication of this value by 
fg ; ^ r .which is found in table 1 to be 2.66, yields ■: 

the -. value' -je-pp ear ing_ in ..column 12,. Column .1:3 is then the . qilm 
of . columns .1.0 . and 12. _ . . . 

"■ 1 •" The vdlues' found in columns 15 and 18 are 'those given ' 
in table -1 corresponding to .a x * p ■*= 1.' Column 19 is the 
value in column 15 less that in- 'column 18.' If Beveral steps 
in the .cUlcUlati ons -are anticipated , (m+i) can be estimated 

to be about 5 so that from table 2 M* ' -is about 58, The 
product '.of ■■ columns 1-4 and 19- is 'listed in Column 20, the ?r •• 
value of M* in column 21. Column 22, the product of a.o‘l-' 
umns 20 and 21, contains one of the terms in equation (37a), 
the equation . from which (/nti ) .. will be calculated. Column 
23 lists'. f x (n). as: taken from table l. r • The -value of ,I r - 

in column -24 was. cal culated-.-ao cording .to the suggestions of 
reference 6 that is, the. effective width ef the sheet was 
assumed .to be .equal to the widths of t^.0 stri/igor. With the 
value, of taken 'from table^l.A the vfeiue. t.o be .filled, in* in 

column 25 can be e.asily computed* • This- is the .second term ■ 
needed for th,e..calo.ul at i on, -..of t -(m+l).: The .third .term :.-i b.-.‘ ob- 

tained by multiplying together . .tho values' found -in .columns 4 
and 13, and is listed in column 26. Subtracting tho value 
listed in- -column 22 -from that .in. column 26, and dividing the 
difference by the v.al.ue in column 26, gives (m+l) 4 , These 
operations' are listed in columns 27" and .28.'.’ It .should bo 
noted how small 'the value in column 22 is as compared to' that 
in column 26. Now (m+l) 3 and (m+l) . are obtained and 
listed in -columns, 29 and 30* .. . 

In order Vo 'calculate £ raaX from equation (34), the' 

values listed in columns 20, 25, and 26 must be multiplied 
by t-he proper function- s; of (m+l). Thus column 31 contains 
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the value of column 25 multiplied by (m+l) , column 32 the 

value of column 26 divided by (m+l) s ,ah'd column 34 the 

value of column 20 multiplied by f(m+l)'. The. value of this 

latter function is taken from figure .7 and is listed In- col- 

umn 33. The sum- of the- values ' in columns 31,, 32, and 34 is 
given in column 35. Column 36 lists ’2w 0 = 2 inches, as as- 
sumed at the "beginning of this numerical example^ column 3T s 
the cross-sectional area of the stringer plus -effective ’■ 
width eomhinati on. This is 0,0524 square inch since' the 
cross-sectional area of the stringer alone is 0,0324 square 
inch, and the thickness of the sheet is 0.010 inch. Column 
38 contains the value of the function f^(s,n) .as taken 

from table 1. Finally, a multiplication and a division, 
carried out in columns 39 and 40 f - yields frEe value of the 
theoretical critical s-train £ max » Isa. the present case 

e max = 0,00245 is obtained .(col. 40), ' ' ^ 

In column 41 the experimental critical strain is 1 i sted”. 
For specimen 31 it is 0. 00190, ,a value definitely smalle,r 
than the one obtained theoretically. On the other hand, 
s = 4 may not be the correct choice. Indeed, figure 10 of 
reference 3 gives n » 3, approximately, if the geometric 
and mechanical properties of specimen 31, and the experimen- 
tal critical strain are used. The choice of s' a 6 should 
be best, therefore, since it yields .n = .20/6.= 3V3 which 
is close to the value 3. . In table 3 th’e entire, calculation 
is repeated' now on the assumption that s = 6, and not 4, 
and all the- other assumptions are kept- unchanged. - ^his time- 
(m+l) becomes, 7*1, fend the theoretical critical strain 
£ max ■'’O. 00217. This is p.bout 14 per.cent higher than the ex- 
perimental strain. The reported-number, (m+l) ' Of ring 
fields involved -In buckling is 9 as given in column 43. 

To make sure that • s = 6 really gives the minimum’ crit- 
ical strain in the inward bulge type instability of this 
monocoque cylinder, the calculations are once more repeated, 
this time on the assumption that s 'S' 8. This assumption 
corresponds to n = 2.5 and gives (m+l)- =-9.54 and 
€ max - 0.00356. .... 


CONCLUSIONS 


In this report the theory of the inward bulge type of 
general instability has been extended to include the effects 
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of the shear strain energy stored in the sheet covering! 

Tables containing numerical- values of the functions involved 
in the. numerical .calculations- ha,ve be:en : prepared which permit 
an easy and- rapid computation of , the maximum strain at fail-: 
ure, : ..Comparison rof -the theory with the experimental results 
obtained at • C.ALC-IT- and P.I-.-B. show an avorag.o deviation ,of : 

10 percent, and a maximum deviation .of 24 percent.. 

The theory as presented in :thi,s report has. the. follow- 
ing shortcomings': * ' *" v 

1. ' The tables 'containing the numerical values of the 

functions involved in the calculation's have beon prepared 
only -for -specimens having the; same- number f 'o£ stringer's atf- 
those- of the GALCIT and P. I.-B,- t est -series. : The tables' can' r 
be easily expanded to include other numb e r s- ' -'of stringers, 
but it is believed that this, work should- be undertaken only , 
after the questions listed under (2) and (3) are solved in a 
satisfactory- manner... • ....... . . : , - 

. . .... , m I 

* . V ' ■ | ** • . . • ’ ’ * i 1 t -- - , 

2. In the development of .the theory ..an .a-seumptipp. was. 

made 'concerning the shearing .rigid.it.y-. . o^ a ejarved. panel of. . 
sheet-- . .wliiph.- ha,s- buc-kled- under, the action of shearing a. tresses. 
This ac sump tio-n should be replaced by data-.-t.o- bo .obtained ; 
from oxporimont, ■ • . ...... #.•«.. .•> . ; <• • 

;3», Because ;of considerations, .of .economy, all- ex-perijaental 
data -avail a.bl e at .present -on general, instability hftve .OJftjMi -y 
obtained ;wit.h spsc-l-m-ens . much. .smaller . t-han _ac,tjial. aircraft 
structures-. The ..scale of-fo.c-t has nocossitat.o.d -the -uao orf- 
.comparatively .weak.- rings, and. .strong qt.r in-ger s ,. « . ;H-orcovor , - 

all test.s have .been conducted with cylinders of approximately 
the same diameter^. ' It- is believed that the range, .of t'esi^d * 
must ho extended, ’ especially v in the direction Of 7 larger di<-' 
ametors, before tho problem of general instability can be 
"considered as finally* -'solved. 

Polytechnic Institute of Brooklyn, r - 

Brooklyn, •’ -K. Y, March '1-944. : “ • 
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APPENDIX I 


Deflected 'Pattern of a 1 Ping * ■ • - - 

• • • •. • * - . • ■■ . *■- ■ ■ ,• '* . *' . 

i 

It can "be shown from the theory of Fourier series that 
axiy ’ symmetrically deflected pattern of a circular ring c'&n- 
he represented by: • 

CO * 

(Al.l) 


. k = i , 2 , 3 , . . . , . j cp j ■ < .Tt'/n_ • 

From 1 the- basic assumption, of ineitensiohal deformations (see 
equation (4gi,"))V it follows that 'Y. 

w t = 6 ^ . .(G^/k) sin(knqx) • . ■ ( A1 . 2 ) 

■ • .. . ■ . k = X 

These functions ■ are ' seen.' tso- satisfy automatically the' five 
conditions, a, b , c , e, ahd f , laid down in equations (4) of 
this paper* " / * • ' • • - . • . . :: . ~ * ~ 

. - . ■■ 

Now "by imposing the requirement, _ 

- •; " (d p w r /dcp :p )(^ F = ; '0"' ' _ (A1.3) 

where .- •?;: - - - ■ -. . ~ 

p = 0, 2, 4, 6,..., 2 (u-1 ) 

. ’ • ' ' •• '• ■ ' - - : . . ' ' ' v , , ‘ . , ’’ : ■ '■ 

the remaining two - requirements (4d) and (-4g)- will "bo ful- 
filled (proyidpd u.j> 2)' and at tho same- timo-j (u~2) addi- 
tional' conditions will arise. \ ' ’ • 

Let 0 1 s* 1 and Oj^. = 0 if k > u + 1. Then this 
criterion. -can be. written as a set of., u . linear equations in 
u -unknowns*. , • . • ' . ....... 


w. 


= - g 


l: 

k=i 


O^oos (kncp) 
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1>* si *f» 1 

1+ o a - o a + ••• + (-1) o k + •••+ (-1) 0. 


u+x 


= 0 


- 1 + 2 3 0 2 - 3 3 0 3 + • • • + (-1) k 2 C v + 


+ [ ( “ ] 


+ | (-1) U+1 (u+i) 3 c u+i = 0 


2 (u-l ) 2(u-i)' . s k a(u-i) 

1+2 ' ’ C s - 3 0 3 + • • • + (-1) • k 


°k + 




1) 


n+ 1 




The common determinant' 'of .this ' syst eni is : 


A — 


2 ° 

2 2 


3 s 


(u+l)° 

(u+l) a 

(u+1) 4 


.gsCu-x), gS(u-l) 


(u + l) a(,1 - V) 


(A1.4) 


(A1 » 5 J 


If u is even, the plus sign is valid' if 11/2 is even, the 
minus if . u/2 . i's .odd’, If u is odd, the "plus 'Sdgn~is 
valid if' (u-l)/2 is ,eVen, the minus if (u-l)/2’ is odd. 

It will now he shown that equation (Al;5) can he writ- 
ten as 'the double ^product J j 


(A1.6) 


u (u+ . ) 

a * * n n (ja _ i a J, j > i 

i=3 j=3 


where the_ capital pi’s signify a douhlo series of multipli- 
cations': 

(u) (u+l ) • • • _ 

n n E(3 B -2 S )][(4 2 -2 S )<4^-3 2 )3[(5 2 -2 2 ) (5 S -3 2 ) (5 3 -4 2 )]* 

i = a 5=3 . • ‘ " ' ■ -■ • • 


[(u+l) S - 2 2 ][(u+l} 2 - 3^3 • »’ C(u+1 ) 2 - u S ] (Al.Sa) 
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This is equivalent to the combination of the difference of 
the squares of all the numbers 8, 3, 4, (u+1) taken two 

at a t ime . : 1 

The (i-l)th column of the determinant consists of the 
consecutive .elements farmed bjr rai sing • the ' number i to all 
even powers from 0 to 2 (u-T) . If this number" -i is re-' 
placed by either of the numbers then the (i-l)th and 

(j-l)th columns will become, id,en ; t'I'C.al« v ’Then two columns of 
a determinant are identical, Its value become s zero. 

Now it is; obvious that when, A ; . ie ,e.xpande£ out by -the 
method of minors, it will take the form:" 



consist s 


where P 
tiers 2, 3,i 
nomial in i 
2(u-l) roots, 
roots of; A..* 
set ’of.' (u-T.), 


only 

‘ i * • 1 

i j +i . 

of' deg- 
But it 


of powers and products of the num- 
*: (u+1).' Consider this as a poly- 
2(u-l). Then it must cqjita.in exactly 
was seen that i = ±3 is a pair of 
■Hence all, the rcjots.must be included in the 
numbers" the values of ~W5ti£clx •• j • ; cah take -on.' ' 

, ■ ■ ■ - - ' ' : < '■ . - : . - f ■ o v ; • r- : • 

'C o r re sp'o riding t o the se' ' rob t s * there wi 1 T ’be 'the factor’s 
(j-i) and . , ( J+i), -or, .what is equivalent to this_set, the 
factors ( i'2) r , "the latter'' U-'l _ in number It is' seen 
that all columns of, a are similarly- const tUc ted in "that v 
they have a constant number raised ,to the same even powerB. 
Thus;all the factors o.f. A will .. lie of the form given by 
(ll . 6) . There s’t i'll' "remaihs- the 'questi’pn "whether (AI.6) 
might be multiplied by some constant factor. If the first 
numbers in' each parenthesis of (A1.6) are multiplied together, 
there results- the term - 


T = (2°)(3 2 )(4 4 )(5 6 )(6 8 ) ... (i 2k ) 


(u+1) 


2_(u-l) 


which is the principal diagonal term of A, which term is 
known to be positive. Consequently (A1.6) represents the 
correct evaluation of. 'A., In an analogous manner, the nu- 
merator determinant for any coefficient, C^, can be put in 
the- form: 


* 


(u) 

AC k - ± n 


i-1 


(u+1 )' 

n 

j=2 


(J 2 -i 2 ). 


j > i. 


J t k 
i k 


(A1.7) 
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where the sign, is the same as - ;that ; of equation (A1^6). Thus 
equation (A1.7) divided "by equation (A1.6) gives 



(u + 1 ) 

/ u+x 




(A1.8) 


c, = n (i a 

’-■i5/ n 

i 2 - k s 

» 

i ^ k 


i = 2 

/ i=a 






In the denominator the absolute values 

were used 

as in- 

dicated by tho bars. 

This is permissible, 

si-nco . Cg 

must 

be positive because. 

. A and Aq^ have 

the 

same sign 

as 

stated. The values 

of the coefficients are 

now calculated 

for 

three cases-: 






For 

u = 2 , C '2 =1,6, - 

C 3 = 0.6 




(A1.9a) 

For 

u = 3, 

0 3 = 9/7, 

C 4 = 2/7 


( A1 . 9b ) 

For 

u = 4 , 03=16/7', 

: 0 3 = 27/14, 

, 0 4 = 16/21, 

C E = 5/42 

( A1 . 9 c ) 

Thus 

the equation 



— 


- - -• 

w r 

= -6 Ceos (nep) + 

1.6 c o s ( 2n<p ) + 0 ; i 

6 cos 

(3ncp)] 

(Al.lOa) 

satisfies all seven 

-conditions 

enumerated in equations (4a) 

to (4f). The function 





w r 

= <-{3 [cos(ncp) + 

(16/7) cos (2ncp) + 

(27/14) cos(3ncp) 


+ (16/31) cos(4nqp) 

+ (6/43) cos 

(5ncp) ] 

(Al.lOb) 


will fulfill two additional requirements at qp = cp 0 . 

Equations (Al.lOa) and (Al.lOb).are plot.ted in figure^2y3 
it is seen how closely the latter curve hugs the original 
oylinder in the region of the termination of the "bulge. it 
is evident thai the more conditions are assumed at the end 
of the "bulge* the smoother the curve will "become. In the 
limit, when an infinite number of derivatives are made to 
vanish at cp = cp 0 , tho deflected shape probably will bocome 

coincident with the original circular ring (except perhaps 
in an interval clbs'e to cp = 0),‘ 

In reality, however,, the deformed pattern will be such 
that only a few of the derivatives vanish at the end point". 

The correct number of terms to be. used should be determined 
from considerations ’ of the minimum of the total potential of 
the system. This cannot be done without the knowledge of tho 
law governing the shearing rigidity of tho buckled panels. 
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APPENDIX II 

CALCULATION OP SUMMATIONS 
Intr oducti on 


In' order to carry out alL sutgeouent summations, 
shall. be made of the fundamental identities: 

cos x + cos(x + a) + cos(x 4* 2a) + ••• + cos [x + (n <- l)a] 

' = cos [x 4* (n - l)a/2] sin(na/2) Csc(a/2) ‘ 


sin x'+ 'sin(x 4- a) + sin_(x + 2a) + * * • + sin [x 4* (n - l)a] 
= sin [x + (n - l)a/2] sin(na/2) csc(a/2) 


PrObf' of these identities' m ky he found in reference 7 
= 0 f> and-. . (n - l) , ig replace.d by N, there result 
orfas' 5 " ‘ •• " • 

N . 

-.? os(la) 7' sin .[(N + l)a/S 3 . osc(.a/2), . • 

" = (l/2') f 0r + cos(-Na) '+ sin(Na) cot(a/2)] ‘ *■ 


i=o . 
N 


7 sia('ia')' 

i=c . . 


s“ SinCNa/2) Sin C(N .'4-1 ) a/2] •• csc(a/2) " ' ■: : 

• s' ( 1 / 2 ) t)l - cos(lla) + sih(Na). cot(a/2)3 
Consequently: : • . * v i • • • •. *. -* n’:. ; * » . 

(s/b)-l ■- ■■■: •• -< "’I : I ■ 

y oos(2TriX/s) = cos[(s/2)'-? l] (tiX/s) 


sin(s/ 2 ) (ttK/s) OscCttX/s) 


i=0 


. =. 0, if X. ia. even.')- 

■. u . • . Vj when. K £ Is* 

- =; 1, if • K . Is. Odd J 


= s/2, 


whpn X.= Is 


use 

(A2.1a) 

(A2.1b> 

.11 

the 

: m 

; (A2.2a) 

(A2.2b) 

i (A2.'3) 
” (A2.3a) 

’ CA2.3o) 
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where K, l, and s are positive integers. Subtracting 1 
from equation (A2.3), 

i 


( 3 / 2 )-! 



y 1 cos(2tt1E/s) = -1, if £ is even"'] 


,(A2.4a) 

L 

= 0, if E is odd J 

where K £ l s 

(A2.4b) 

= (s/2)-l, 

when E = Is 

(A2.4c) 


Also from (A2,la)j 


( s/2 )— 1 

cos(2ttK/s) (i + l/2) = coe [(ttK/s) 

i=o 

+ (s/2 ~ l) (ttK/s)] sin (s/2) (ttK/s) csc(ttE/s) 

= 0 for all K, K / Is (A2.5a) 

= (-1) (s/2) for K = Is .. (A2.5'b) 


Terms Involved in the Strain Energy 
Stored in the Stringers 

With the aid of equations (A2,4a) to (A2.4c) of this 
appendix, the summations appearing in equation (l?) of this 
report can now be calculated. 
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bo 


(s/a )- 1 

256 ^str r + 2 I s tr r [2 coe(Tri/s) + 9 oos(jni/s )+5 cos(5tt1/s) ] 

i=i 

(s/a)-i 

+ (2/n 2 ) I ^ [4 sin(TTi/s) + 6 sinCjrTi/s) + 2 sin.(5rTi/s ] s 

str t / . 


i=i 


(s/a)-i 


— 256 Igtr r + ^ ^str r 


+ 21 


C2 3 oos a (TTi/s) + 9 2 cos 2 (3ttI/s) 

i=i 

+ 5 2 CO 5 ( 5 ^ 1 / s)] 

str r ^ 2 t(2)(9) cos (ttI/s) cos(3iTi/ , s) 

+ (2)(5) cos{ni/s) cos(5fTi/s) + (9)(5) cos(3ni/s) cos(5rri) 

+ (2/n 2 )l str t [4 2 sin 2 (Tri/s) + 6 2 sin. a (3ni/s) + 2 2 sin. 2 (5TTi/s) ] 

+ (2/n a )l str t ^ 2 [(4) (6) sinCfTi/s) sin(3TTi/s) 

+ (4) (2) sinOrri/s) sin.(5Tri/s) + (6)(2) sia(3Tri/s) sin(5TTi/s) ] 


s I 


str r 


<{ 256+2 £ (l/2)(2 2 + 9 Z +5 S ) + 2 ^ (1/2) C2 a cos(2iTi/s) 

r 

+ 9 2 cos(6ni/s) + 5 2 cos 10rri/s3 

2 ^ C(2)(9)(cos 4ni/s+ cos 2-rri/s) + (2)(5)(cos 6rri/s + cos4ni/s) 

< + (9) (5)(cos 8rri/s + cos2TTi/s)3^ 

(Istr t/n 2 ) {2 £ d/2)(4 2 + 6 2 + 2 3 ) - 2 £ (l/2) [4 2 cos(2Tri/s) 

+ 6 2 cos(6tt1/s) + 2 s cos(lOrci/s) 3 
2 £ C(4)(6)( cos 2ni/s - cos 4rri/s) + (4) (2) (cos 4tt1/s - cos 6rri/s) 
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All these summations fall directly under the types of equa 
tions (A2.4a) to (A2.4c). By use of these formulas, there 
is obtained: 


* rtr C.,n) = I otr r {t(«/2) - 1 ) fa 3 + 9 2 + 5 8 ) - 2 C(2)(9) + Oi <5> 

+ (9) (5)3 + 256 

+ (i, tr t /u 2 ) -f" <(s/2) - >x ) (4. 2 + 6 2 +2 2 ) t2[ft)(6)-W(2)t(6)(2)]} 
if s > 4. Or 

^str ( s » n ) = (55s) I s tr r (2Ss/n )lgtr t • 

= fgtr r( s > EL ) ^str r "t ^str t^ s >^) ^str t 


htM.' ^ 


^ * 7 + 


(A2.6a) 


It follows from actual calculations (or equation (A2.4c) ) 
that 


F s tr(s.n) = 400 I sfcr r + (64/n 2 )l str t if s = 4 (AS. 6b) 


and 


•® I str( s » 11 ) — 256 Igtr r 


if s = 2 (A2.6 c) 


Terms Involved in 
St or ed 

The summation 'over i 
to be evaluated. 


the Shear Strain Energy 
in the Sheet 

appearing in equation (23) has 
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(«/a)-i 

^ Csln(rd/g)+ sln.TT (t+ l)/s3 + 6 C sin(3ni/s) + sin 3 tt(1+1)/s 3 
i-o a 

+ 2 £Bin( 5 rTi/s) + ein 5^(1 + l)/s3j- 

s Z 16 { (i/,2) El - cos(2tt1/s)3 + Ccos(tt/s) - eos(2rr/s)(i + 1/2)3 

+ ( 1/2 Cl - oos( 2 ir/s)(i + 1 )]^ 

+ 36 ^(l/ 2 ) [l — cos(&ni/s )3 + CcosCjrr/g) - cos( 6 tt/s )(1 + l/ 2 )} 

+ (l/ 2 ) [l - cos(Gn/s)(i + l)]^ 

+ V 4 £(l/2) [1 - COB (lOrrl/s) ] + £ 008 ( 5 ^/*) - cob(10tt/s)(1 + l/2)3 

+ (l/2) Cl - co«(lOn/s)(i + l) 

+ ^ (4) (6) -^Ccos^i/s) - oos(4rA/s)3 + C cob 2 tt (1 - l/2)/s 

- cob(4it/s)( 1>+ l/4)] + Ccoa(2iT/s)(i + 3/2) - ooB(4rc/a)(i + 3/4)3 

+ Coob(2tt/s)(1 + 1) - oos(4tt/s)( 1 + l)3^| 

+ ^ (4) ( 2 ) ^£co8(4ni/s) - co8(Gni/s)3 + CcoB(4n/a)(l + 5/4-) 

- C08(&n/s)(i 4- 5/6)] + £co6(4tt/s)(1 “ 1/4-) - ccs(&rr/s)(l + l/6)3 

+ £ cos(4rr/8) (i + l) - cob( 6 tt/s )(1 + l) ]j 
+ ( 6 ) (2) ^Co< 5 b( 2 tt 1 / 8 ) - cos(8rTi/s)3 + £ oos(2tt/s) (i + 5/2) 

- co8(Stt/s)(1 + 5/^)3 + Coos(2n/s)(l - 3 / 2 ) - oos( 6 n/s)(l + 3/8)3 

+ Ceo8(2TT/8)(i +- l) - cos(Sn/s)(i + l)3^ 

a £ |U6 + 36 + 4 ) + £l6 cos(tt/s) + 36 cos(3rr/s) + 4 cos( 5 rc/s )3 

+ £fi(«) co«(2ti/«) ( 1+1/2) + f*(«) 008 (Wb)( 1 + 1 / 2 ) +«*(«) cot(6ir/«)(i+l/2) 
+ f 4 (e) cos(8TT/«)(i + 1/2) + ^(«) oot(l0n/8)(i + 1/2)3 
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where the f’s are functions of s alone. 

from equation (A2,5a) it follows immediately that this 
summation reduces to 

) cos 3tt/s + cos 5rr/s) J s > 4 (A2.7a) 

[ 56 +4 (4 costt/s 

+ 9 cos Jw / s + cos ^/s)] s > 4 (A2.7b) 



Carrying out the summation numerically, it is found that 


a o f ish (s,n) = ^(128/^) 


if s = 4 (A2.7c) 


= 0 


if s * 2 


(A2.7d) 


There is left the evaluation of the last bracket in 
•equations (23). 

'(m+i)-i 



[cos 2-rr(j + l)/(m + l) - cos(2rr/m + l)j ] 2 
{(1/2) [1 - cos 4rr(j + 1 )/ (m+ l) ] + (l/2) [l- cos 4 ttJ (m+ l) ] 


- cos(4rr/m+ l) (j + l/2) 


cos (2Tr)/(m+ l) J- 


s (m + l) [l - cos(2Tr)/(m + l)3 - cos 4rT/(m.+ l) 

(m+i)- 1 

7 [cos(4rr/(m + l)] [ j + (l/2)] 

J 

j=0 

from equation (A2.5a) it follows that the last bracket van- 
ishes when (m + l) > 2, Hence the summation reduces to: 

f (m+ l)/(m+ l) = (m+l) [l-cos 2n-/(m+l)3 for (m+ l) > 2 (A2.Sa) 

f(m+l)/(m+l) =16 for (m+l) = 2 (A2.Sb) 




NACA TN NO. 968 


Table 1. Strain Energy Expressions 




— ■ 1 1 I I ■ ■ — — 

Mii«a 



s 

n 

f str.r( s » n ) 


fiahCsjn) 

f 2ah(“» n ) 

p - 1 

*2ah( s » n ) 

p - 2 

f r ( n ) 

f^ejn) 

10 

2.5 

400 

10.24 

s i> 4 
20.480 

15.758 

17.837 

20016.4 

380.783 

12 

3 

400 

7 1/9 

‘ 14.222 

11.897 

12.969 

35445 • 3 

386.866 

16 

■ 

400 

4 

8.000 

7.249 

7.609 

86071.5 

392.734 


5 

400 . 

2.56 

5.120 

4.809 

4.961 

169986.4 

395.387 



330 

60.48 

■21 

MH 



■ 

m 


12 

<P 

330 

42 

99.590 

80.229 

88.993 

9434 

333.765 

16 

8/3 

330 

23.625 

56.020 

49.745 

52.728 

24529.8 

332.915 

20 

10/3 

330 

15.12 

35.852 

33.255 

34.513 


332.088 

24 

■ 

330 

10.5 

24.898 

23.637 

24.254 

86071.5 


16 

2 

440 

56 

s ■ 8 
166.055 

138.426 

151.153 

9434 

433-576 

20 

2.5 

440 

35.84 

106.275 

93.259 

100.248 

gg 

436.819 

24 

3 

440 

24 8/9 

73.802 

68.204 

70.913 


BE 

*. 

40 

5 

440 

8.96 

26.569 

25.834 

26.197 

169986.4 

439.527 




































































































20 2 


24 


f str.r( s » n ) 

f str.t( s » n ) 

550 

70 

550 

48.61 

550 

- 17.5 



195.003 

211.672 

71.409 ? 

75-587? 

55-415 

. 56.588 



Table la 


2 3 4- 5 6 7 8 16 


| 39.478 9.869 4.386 2.467 1.579 1.097 0.8057 0.6169 0.1542 



54.9 58.4 60.3 61.6 62.3 62.9 63.2 64-9 
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Table 3- 


Calculation of the Crit:LcalCompressive Strains 


Cl] 
GAicrr 
Spec . No . 


[2] [3] C43 , 

S L x (2^1-1)' 

(iK.) (in.-2) 


40 1 e 


40 16' 

20 8 


10 8 


Mil 


50 

10 

51 

10 

5-1 

10 

54 

40 


K»1 


Poly. 

Specs. 


Thick 

Sheet 


9.870 

2.467 




[63 

2w« 

(in.) 


[73 [8] [93 

r t Istr. 
(in.) (in.) (in4xl 


15.92 

15.92 

15.92 

15.92 

15.92 


a^) 


15.92 

15.92 

15.92 

15.92 


\Wt3EEM 

tmwiAi 


1.5 15.92 0. 

1.5 15.92 0. 

1.5 15.92 0.015 

2 15.92 0.015 


2 

2 

2 

2 

2 


1 


15.92 

15.92 

15.92 

16.0 




iir.i 


nmMi 


16.0 

0.010 

16.0 

0.010 

16.0 

0.010 

16.0 

0.010 


16 

5 

1.579 

6 

16 

5 ' 

1.579 

6 


[103 

(9) 

x ^str.r( s > n ) 


RSI 


2520 

.831 

2190 

.723 
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XVr** * f {tp* » •teo'jxyf.iCii. 
- 

Table 3. Qal relation of the Critical Compressive Strains 



[i] 

[11] 

5S 

X fstr.t(s,n) 

[13] v 

[14] 

[15] 

[16] 

[17] . 

GALcrr 
Spec. No. 

Iatist , 
(intafclO -6 ) 

Fstr( s J n ) 0 

(10)+ (12) 

.0307tr/S 

OnfxlO-6) 

fl.sh(s,n) 

a l 


25 

337 

.0590 

0.407 

122 


0.90 

■ 

26 

337 

.0590 

0.407 

122 

Irivln 

0.90 


27 

337 

0.102 

0.519 

122 


0.95 


28 

723 

0.065 

0.359 .v» 
0 . 329 ■*' 

122 


0.80 

■ 9B9| 

30 

723 

0.050 8< 

244 


■H 


SPH 

■PA«i:aaai 

.01093^' 

0.277? 

244 

HWUSI 


1 

kS 

408 

0.0616^'*/ 

0.276 

244 



1 

H 

408 

0.1465 ,S ^ 1 ’ 

0.432 

244 

llKaBSli 


r 

^ 32 

408 

0.0616 

0.276 

244 

35.85 

1 

1 

35 ' 

942 

.0965 

.369 

489 

20.48 

1 

2 

36 

~“W 

75955 

39 

489 

20.48 

1 

2 

37 

942 

' .0965 

.369 

489 

20.48 

1 

2 

38 

408 

.0616 

0.276 

244 

35.85 

1 

1 

3^ 

408 * 

.0616 

0.276 

244 

35.85 

1 

2 


99 

.0300 

0.391 

122 

105-78 

0-95 

1 

39 

— 99 

75305 

0.391 

122 

105.78 

0.95 

1 

41 

478 ' 

.0836 

0.480 

183 

57.74 

0.875 

1 

42 

478 

.0836 

0.480 

183 

57-74 

0,90 

1 

43 

478 

• 144 

0.619 

183 

57.74 

0.90 

1 

45 

1117 

.368 

0.419 

366 

35.85 

1 

1 

45 

TH7 

.168 

.0.419 

366 

35.85 

1 

2 

46 

1117 

.168 

0.419 

• 366 

35-85 

1 

1 

„ 46 

1117 

.168 

0.419 

366 

35.85 

1 

2 

47 

1117 

.168 

0.419 

366 

35-85 

1 

1 

47 

1117 

.168 

0.419 

366 

35.85 

1 

2 

49 

“2515 

7355 

0.5l5 

733 

20.48 

1 

2 

50 

2010 

.206 

0.515 

733 

20.48 

1 

2 

51 ■ 

2010 

.206 

0.515 

733 

20.48 

1 

1 

51 

2010 

.206 

0.515 

733 

20.48 

1 

2 

54 

140 

.0673 

0.510 

123 

174 

1 

1 

55 

58 

U5 

337 

— OTTJ 
oTisr 

0.606 

0.399 

123 

246 

265-9 

106.3 

1 

1 

1 

1 

58 

337 

0.121 

0.399 

246 

106.3 

1 

2 

59 

337 

0.121 

0.399 

246 

106.3 

1 

1 

59 

337 

0.121 

0.399 

246 

106.3 

1 

2 . 

61 

337 . 

0.121 

0.399 

246 

15573 

1 

1 

61 

337 

0.121 

0.399 

246 

106.3 

1 

2 

63 

140 

.0673 

0.510 

123 

174.0 

1 

1 

64 

140 

.100 

0.606 

123 

265.9 

1 

1 

'65 

140 

.0348 

0.287 

128 

73-80 

1 

1 

-■66 

140 

.0348 

0.287 

128 

73.80 

1 

1 

67. 

408 

.0290 

0.283 

256 

14.22 

1 

1 

68 

337 

.0241 

0.273 

256 

14.22 

1 

1 

Poly.Spe 

cs. 







Thick 

Sheet 

728 

.172 

1.00 

384 

56.02 

1 

1 

Thin v 
Sheet 

509 

.120 

.843 

■384 

56.02 

1 

1 
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FIG. Id DISTORTIONS OF A RING AT FIG. lb. COMF^RISON OF DEFLECTED 

BUCKLING. PATTERNS OF A RING. 

FORMER ASSUMED SHAPE. 

REVISED SHAPE. 

FOR n ~ 3 
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A 


B* C' 



FIG. 4. SHEAR DISTORTION OF PANEL. 





FIG. 5. EFFECTIVE SHEAR RIGIDITY VARIATION. 

g'-g (i- cos 40 

g 1 - g C»— cos ♦) 

* — g'« g (i -Jcos $ ) 
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Flga. 6,7 




FIG. 7. FUNCTION OF NUMBER OF RING FIELDS INVOLVED IN FAILURE. 


